Author Correction: Multiplexed droplet single-cell RNA-sequencing using natural genetic variation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Functional CRISPR dissection of gene networks controlling human regulatory T cell identity.

Human regulatory T (Treg) cells are essential for immune homeostasis. The transcription factor FOXP3 maintains Treg cell identity, yet the complete set of key transcription factors that control Treg cell gene expression remains unknown. Here, we used pooled and arrayed Cas9 ribonucleoprotein screens to identify transcription factors that regulate critical proteins in primary human Treg cells under basal and proinflammatory conditions. We then generated 54,424 single-cell transcriptomes from Treg cells subjected to genetic perturbations and cytokine stimulation, which revealed distinct gene networks individually regulated by FOXP3 and PRDM1, in addition to a network coregulated by FOXO1 and IRF4. We also discovered that HIVEP2, to our knowledge not previously implicated in Treg cell function, coregulates another gene network with SATB1 and is important for Treg cell-mediated immunosuppression. By integrating CRISPR screens and single-cell RNA-sequencing profiling, we have uncovered transcriptional regulators and downstream gene networks in human Treg cells that could be targeted for immunotherapies.

Genetic analysis of isoform usage in the human anti-viral response reveals influenza-specific regulation of ERAP2 transcripts under balancing selection.

While genetic variants are known to be associated with overall gene abundance in stimulated immune cells, less is known about their effects on alternative isoform usage. By analyzing RNA-seq profiles of monocyte-derived dendritic cells from 243 individuals, we uncovered thousands of unannotated isoforms synthesized in response to influenza infection and type 1 interferon stimulation. We identified more than a thousand quantitative trait loci (QTLs) associated with alternate isoform usage (isoQTLs), many of which are independent of expression QTLs (eQTLs) for the same gene. Compared with eQTLs, isoQTLs are enriched for splice sites and untranslated regions, but depleted of sequences upstream of annotated transcription start sites. Both eQTLs and isoQTLs explain a significant proportion of the disease heritability attributed to common genetic variants. At the ERAP2 locus, we shed light on the function of the gene and how two frequent, highly differentiated haplotypes with intermediate frequencies could be maintained by balancing selection. At baseline and following type 1 interferon stimulation, the major haplotype is associated with low ERAP2 expression caused by nonsense-mediated decay, while the minor haplotype, known to increase Crohn's disease risk, is associated with high ERAP2 expression. In response to influenza infection, we found two uncharacterized isoforms expressed from the major haplotype, likely the result of multiple perfectly linked variants affecting the transcription and splicing at the locus. Thus, genetic variants at a single locus could modulate independent gene regulatory processes in innate immune responses and, in the case of ERAP2, may confer a historical fitness advantage in response to virus.

Genetic determinants of co-accessible chromatin regions in activated T cells across humans.

Over 90% of genetic variants associated with complex human traits map to non-coding regions, but little is understood about how they modulate gene regulation in health and disease. One possible mechanism is that genetic variants affect the activity of one or more cis-regulatory elements leading to gene expression variation in specific cell types. To identify such cases, we analyzed ATAC-seq and RNA-seq profiles from stimulated primary CD4+ T cells in up to 105 healthy donors. We found that regions of accessible chromatin (ATAC-peaks) are co-accessible at kilobase and megabase resolution, consistent with the three-dimensional chromatin organization measured by in situ Hi-C in T cells. Fifteen percent of genetic variants located within ATAC-peaks affected the accessibility of the corresponding peak (local-ATAC-QTLs). Local-ATAC-QTLs have the largest effects on co-accessible peaks, are associated with gene expression and are enriched for autoimmune disease variants. Our results provide insights into how natural genetic variants modulate cis-regulatory elements, in isolation or in concert, to influence gene expression.

Author Correction: Discovery of stimulation-responsive immune enhancers with CRISPR activation.

In this Letter, analysis of steady-state regulatory T (Treg) cell percentages from Il2ra enhancer deletion (EDEL) and wild-type (WT) mice revealed no differences between them (Extended Data Fig. 9d). This analysis included two mice whose genotypes were incorrectly assigned. Even after correction of the genotypes, no significant differences in Treg cell percentages were seen when data across experimental cohorts were averaged (as was done in Extended Data Fig. 9d). However, if we normalize the corrected data to account for variation among experimental cohorts, a subtle decrease in EDEL Treg cell percentages is revealed and, using the corrected and normalized data, we have redrawn Extended Data Fig. 9d in Supplementary Fig. 1. The Supplementary Information to this Amendment contains the corrected and reanalysed Extended Data Fig. 9d. The sentence "This enhancer deletion (EDEL) strain also had no obvious T cell phenotypes at steady state (Extended Data Fig. 9)." should read: "This enhancer deletion (EDEL) strain had a small decrease in the percentage of Treg cells (Extended Data Fig. 9).". This error does not affect any of the main figures in the Letter or the data from mice with the human autoimmune-associated single nucleotide polymorphism (SNP) knocked in or with a 12-base-pair deletion at the site (12DEL). In addition, we stated in the Methods that we observed consistent immunophenotypes of EDEL mice across three founders, but in fact, we observed consistent phenotypes in mice from two founders. This does not change any of our conclusions and the original Letter has not been corrected.

A Calmodulin C-Lobe Ca2+-Dependent Switch Governs Kv7 Channel Function.

Kv7 (KCNQ) voltage-gated potassium channels control excitability in the brain, heart, and ear. Calmodulin (CaM) is crucial for Kv7 function, but how this calcium sensor affects activity has remained unclear. Here, we present X-ray crystallographic analysis of CaM:Kv7.4 and CaM:Kv7.5 AB domain complexes that reveal an Apo/CaM clamp conformation and calcium binding preferences. These structures, combined with small-angle X-ray scattering, biochemical, and functional studies, establish a regulatory mechanism for Kv7 CaM modulation based on a common architecture in which a CaM C-lobe calcium-dependent switch releases a shared Apo/CaM clamp conformation. This C-lobe switch inhibits voltage-dependent activation of Kv7.4 and Kv7.5 but facilitates Kv7.1, demonstrating that mechanism is shared by Kv7 isoforms despite the different directions of CaM modulation. Our findings provide a unified framework for understanding how CaM controls different Kv7 isoforms and highlight the role of membrane proximal domains for controlling voltage-gated channel function. VIDEO ABSTRACT.

Multiplexed droplet single-cell RNA-sequencing using natural genetic variation.

Droplet single-cell RNA-sequencing (dscRNA-seq) has enabled rapid, massively parallel profiling of transcriptomes. However, assessing differential expression across multiple individuals has been hampered by inefficient sample processing and technical batch effects. Here we describe a computational tool, demuxlet, that harnesses natural genetic variation to determine the sample identity of each droplet containing a single cell (singlet) and detect droplets containing two cells (doublets). These capabilities enable multiplexed dscRNA-seq experiments in which cells from unrelated individuals are pooled and captured at higher throughput than in standard workflows. Using simulated data, we show that 50 single-nucleotide polymorphisms (SNPs) per cell are sufficient to assign 97% of singlets and identify 92% of doublets in pools of up to 64 individuals. Given genotyping data for each of eight pooled samples, demuxlet correctly recovers the sample identity of >99% of singlets and identifies doublets at rates consistent with previous estimates. We apply demuxlet to assess cell-type-specific changes in gene expression in 8 pooled lupus patient samples treated with interferon (IFN)-ß and perform eQTL analysis on 23 pooled samples.

Discovery of stimulation-responsive immune enhancers with CRISPR activation.

The majority of genetic variants associated with common human diseases map to enhancers, non-coding elements that shape cell-type-specific transcriptional programs and responses to extracellular cues. Systematic mapping of functional enhancers and their biological contexts is required to understand the mechanisms by which variation in non-coding genetic sequences contributes to disease. Functional enhancers can be mapped by genomic sequence disruption, but this approach is limited to the subset of enhancers that are necessary in the particular cellular context being studied. We hypothesized that recruitment of a strong transcriptional activator to an enhancer would be sufficient to drive target gene expression, even if that enhancer was not currently active in the assayed cells. Here we describe a discovery platform that can identify stimulus-responsive enhancers for a target gene independent of stimulus exposure. We used tiled CRISPR activation (CRISPRa) to synthetically recruit a transcriptional activator to sites across large genomic regions (more than 100 kilobases) surrounding two key autoimmunity risk loci, CD69 and IL2RA. We identified several CRISPRa-responsive elements with chromatin features of stimulus-responsive enhancers, including an IL2RA enhancer that harbours an autoimmunity risk variant. Using engineered mouse models, we found that sequence perturbation of the disease-associated Il2ra enhancer did not entirely block Il2ra expression, but rather delayed the timing of gene activation in response to specific extracellular signals. Enhancer deletion skewed polarization of naive T cells towards a pro-inflammatory T helper (TH17) cell state and away from a regulatory T cell state. This integrated approach identifies functional enhancers and reveals how non-coding variation associated with human immune dysfunction alters context-specific gene programs.

CRISPR/Cas9-mediated PD-1 disruption enhances anti-tumor efficacy of human chimeric antigen receptor T cells.

Immunotherapies with chimeric antigen receptor (CAR) T cells and checkpoint inhibitors (including antibodies that antagonize programmed cell death protein 1 [PD-1]) have both opened new avenues for cancer treatment, but the clinical potential of combined disruption of inhibitory checkpoints and CAR T cell therapy remains incompletely explored. Here we show that programmed death ligand 1 (PD-L1) expression on tumor cells can render human CAR T cells (anti-CD19 4-1BBζ) hypo-functional, resulting in impaired tumor clearance in a sub-cutaneous xenograft model. To overcome this suppressed anti-tumor response, we developed a protocol for combined Cas9 ribonucleoprotein (Cas9 RNP)-mediated gene editing and lentiviral transduction to generate PD-1 deficient anti-CD19 CAR T cells. Pdcd1 (PD-1) disruption augmented CAR T cell mediated killing of tumor cells in vitro and enhanced clearance of PD-L1+ tumor xenografts in vivo. This study demonstrates improved therapeutic efficacy of Cas9-edited CAR T cells and highlights the potential of precision genome engineering to enhance next-generation cell therapies.

A tissue checkpoint regulates type 2 immunity.

Group 2 innate lymphoid cells (ILC2s) and CD4+ type 2 helper T cells (TH2 cells) are defined by their similar effector cytokines, which together mediate the features of allergic immunity. We found that tissue ILC2s and TH2 cells differentiated independently but shared overlapping effector function programs that were mediated by exposure to the tissue-derived cytokines interleukin 25 (IL-25), IL-33 and thymic stromal lymphopoietin (TSLP). Loss of these three tissue signals did not affect lymph node priming, but abrogated the terminal differentiation of effector TH2 cells and adaptive lung inflammation in a T cell-intrinsic manner. Our findings suggest a mechanism by which diverse perturbations can activate type 2 immunity and reveal a shared local-tissue-elicited checkpoint that can be exploited to control both innate and adaptive allergic inflammation.

Determination of ubiquitin fitness landscapes under different chemical stresses in a classroom setting.

Ubiquitin is essential for eukaryotic life and varies in only 3 amino acid positions between yeast and humans. However, recent deep sequencing studies indicate that ubiquitin is highly tolerant to single mutations. We hypothesized that this tolerance would be reduced by chemically induced physiologic perturbations. To test this hypothesis, a class of first year UCSF graduate students employed deep mutational scanning to determine the fitness landscape of all possible single residue mutations in the presence of five different small molecule perturbations. These perturbations uncover 'shared sensitized positions' localized to areas around the hydrophobic patch and the C-terminus. In addition, we identified perturbation specific effects such as a sensitization of His68 in HU and a tolerance to mutation at Lys63 in DTT. Our data show how chemical stresses can reduce buffering effects in the ubiquitin proteasome system. Finally, this study demonstrates the potential of lab-based interdisciplinary graduate curriculum.

Generation of knock-in primary human T cells using Cas9 ribonucleoproteins.

T-cell genome engineering holds great promise for cell-based therapies for cancer, HIV, primary immune deficiencies, and autoimmune diseases, but genetic manipulation of human T cells has been challenging. Improved tools are needed to efficiently "knock out" genes and "knock in" targeted genome modifications to modulate T-cell function and correct disease-associated mutations. CRISPR/Cas9 technology is facilitating genome engineering in many cell types, but in human T cells its efficiency has been limited and it has not yet proven useful for targeted nucleotide replacements. Here we report efficient genome engineering in human CD4(+) T cells using Cas9:single-guide RNA ribonucleoproteins (Cas9 RNPs). Cas9 RNPs allowed ablation of CXCR4, a coreceptor for HIV entry. Cas9 RNP electroporation caused up to ∼40% of cells to lose high-level cell-surface expression of CXCR4, and edited cells could be enriched by sorting based on low CXCR4 expression. Importantly, Cas9 RNPs paired with homology-directed repair template oligonucleotides generated a high frequency of targeted genome modifications in primary T cells. Targeted nucleotide replacement was achieved in CXCR4 and PD-1 (PDCD1), a regulator of T-cell exhaustion that is a validated target for tumor immunotherapy. Deep sequencing of a target site confirmed that Cas9 RNPs generated knock-in genome modifications with up to ∼20% efficiency, which accounted for up to approximately one-third of total editing events. These results establish Cas9 RNP technology for diverse experimental and therapeutic genome engineering applications in primary human T cells.

Structural and functional effects of Cu metalloprotein-driven silver nanoparticle dissolution.

Interactions of a model Cu-metalloprotein, azurin, with 10-100 nm silver nanoparticles (NPs) were examined to elucidate the role of oxidative dissolution and protein interaction on the biological reactivity of NPs. Although minimal protein and NP structural changes were observed upon interaction, displacement of Cu(II) and formation of Ag(I) azurin species under aerobic conditions implicates Cu(II) azurin as a catalyst of NP oxidative dissolution. Consistent with NP oxidation potentials, largest concentrations of Ag(I) azurin species were recorded in reaction with 10 nm NPs (>50%). Apo-protein was also observed under anaerobic reaction with NPs of all sizes and upon aerobic reaction with larger NPs (>20 nm), where NP oxidation is slowed. Cu(II) azurin displacement upon reaction with NPs was significantly greater than when reacted with Ag(I)(aq) alone. Regardless of NP size, dialysis experiments show minimal reactivity between azurin and the Ag(I)(aq) species formed as a result of NP oxidative dissolution, indicating Cu displacement from azurin occurs at the NP surface. Mechanisms of azurin-silver NP interaction are proposed. Results demonstrate that NP interactions not only impact protein structure and function, but also NP reactivity, with implications for targeting, uptake, and cytotoxicity.